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ABSTRACT 

In this paper, we present the anti-synchronization of two chaos -exhibiting systems-Bullard and Rikitake 
dynamos using nonlinear active control techniques in a master-slave topology. Nonlinear active control laws 
were derived and added to the algebraic structure of the Bullard slave system and the Lyapunov stability 
criteria was applied to verify the negative definiteness of the error dynamics as a condition for 
antisynchronization of the two systems. Simulation results confirmed the effectiveness of the approach in 
coupling the dynamics of the systems. 
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I. INTRODUCTION 

During the past half of a century since the discovery of chaos in weather by the physicist, Edward 
Lorenz [1] and the experimental control of chaos by other physicists such as Ott et al [2] and Pyragas [3], 
intense research has continued to spurn out literature on the ubiquity of chaos in various natural and man- 
made systems such as economics [4], psychology [5], ecology [6], food [7], . Applications of chaos real life 
designs have continued to grow even as new chaotic systems have continued to be evolved through research 
over the years. Chaos is a phenomenon of nonlinear systems which are extremely sensitive to perturbation in 
their parametric structures leading to unpredictability in their future evolution. Chaos synchronization or anti- 
synchronization occurs when two dissipative chaotic systems are coupled such that, in spite of the exponential 
divergence of their nearby trajectories, synchrony or anti-synchrony is achieved in their chaotic behaviours 
as t — >■ oo . Different methods have been developed to synchronize or anti-synchronize chaotic systems. These 
include Fourier series expansion [8], adaptive control [9], active control [10], fuzzy control [11], impulsive 
control [12], and sliding mode control [13] among others. A simplified architecture for anti-synchronization of 
two chaotic systems is depicted in Fig. 1. The nonlinear control laws synchronize or anti-synchronize the time 
series trajectories of the two coupled systems according to some laws governing the choice synchronization 

scheme where x(0) = x l (0), x 2 (0), x 3 (0) * y(0) = .V[(0), ;y 2 (0), ,v 3 (0) . Synchronization has found 

application in secure communication and other systems while studies on an ti -synchronization of chaos have 
been found to be useful in mitigating electrical power outage and system management [14]-[15]. 
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Figure 1 : Simplified architecture for antisynchronization of two chaotic systems 
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D. DESCRIPTION OF THE RIKITAKE AND BULLARD DYNAMO SYSTEMS 

The Rikitake system is is a simplified dynamic model which attempts to explain the irregular polarity 
switching of the earth's geomagnetic field. The physics of the Rikitake system has been studies by various 
authors [16]-[17], while other authors have studies synchronization of the system using various methods [18] - 
[19]. Mathematically, the system is represented by three coupled-differential equations given as follows: 

x = x.x, - n x , - \jf x 2 
x } = 1 - x { x 2 

(1) 

Where x , x , x } are state variables, y/ , n > 0 are positive constants. For values of y/ = 2, ;r = 3.4641 
the system evolves the 2-D phase portraits and state trajectories shown in Fig. 2 and Fig. 3.. 
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Figure 2. 2D Phase portraits of the Rikitake dynamo system 
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Figure 3. State trajectories of the Rikitake dynamo system 

The Bullard dynamo system is also related to the dynamic problems associated with the earth's 
core and was first studied in detail by E.C. Bullard [20]. The physics of the system may be found in [21]-[22]. 
The algebraic structure of the system is represented by three coupled differential equations given as: 

y 2 = £ - r y 2 - <py] 
y } = ®y t - ry 3 

(2) 
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Where y 1 , y , y are state variables, ® , £ , y , <p > 0 are positive constants. For values of 

0 = 2, £ = 20, y = 2, <p = 4 , the system evolves the 2-D phase portraits and state trajectories shown in 

Fig. 4 and Fig. 5. 
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Figure 4: 2D Phase portraits of the Bullard dynamo system 




0 20 40 SO 0 20 40 60 0 20 40 60 



t(s) t(s) t(s) 

Figure 5: State trajectories of the Bullard dynamo system 

III. ANTISYNCHRONIZATION OF THE RIKITAKE AND BULLARD SYSTEMS 

In this section, the anti-synchronization problem is presented. While the design of nonlinear 
controllers for antisynchronization of chaotic system is not a new challenge, it is however interesting as it has 
continued to elicit interest from researchers due to its potential applications in engineering and non- 
engineering systems. It can therefore be safely said that every (anti)synchronizable chaotic system is a potential 
candidate for multi-dimensional applications in such fields as biomedical engineering, astronomy, economics 
and finance, automotive engineering, radar engineering amongst numerous other fields. Let eq. (1) be the 
master system and eq. (2) be the slave system. The controlled slave system can be represented as 

?i = y^yi - !) - ®y 3 + "i 

y 2 = £ - y y 2 - <p y] + u 2 

y } = &y, - r y 3 + « 3 

(3) 

Where " 1 ' U 2 ' u 3 are nonlinear controllers to be designed. Let the antisynchronization error be represented as 

e = (e t , e 2 , ...e n ) = y + x 
(4) 
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Where x = {x l , x 2 ,...x n ) T , y = (y 1 ,y 2 ,--y„ V and 

e = y + x ; e 2 = y 2 + x ; e 3 = y + x 3 \ e = y n + x . The objective of the study is to design the 
nonlinear active controllers u., i = 1,2,3 such that the trajectories y , y , y of the slave system can come 
into anti-synchrony with trajectories x lt x,, x 3 of the master system, subject to different initial conditions 
x(0) * y(0) such that 

lim ||e|| = lim \\y(t, y (0)) + x(t, *(0)) II = 0 

(5) 

By adding (3) to (1) using the relation in (4), the error dynamics becomes 

e, = y,y 2 - y, - ©y 3 + x 2 x 3 -yr x { + «, 

e 2 = ^- yy,-(py 2 l + x t x 3 - n x, - y/ x 2 + u 2 

(6) 

e 3 = ©y, - yy 3 + 1 - x t x 2 + u 3 
By simplifying (6) and using (4), the error dynamics becomes 

e i = y [ y 2 ~ e i~® e i + x 2 x ?, + ^ + ®~y / } x i + u i 

e 2 = f - ye 2 + yx 2 - q>y\ + x l x i - n e, + n y i - y/ e 2 + y/ y 2 + u 2 

(7) 

e 3 = © e l — © x l — ye 3 + y x 3 + 1 - x^x 2 + u 3 



From (7), the nonlinear active control laws are given as 

"i = - y l y 2 + e l + @e 3 -x 2 x 3 -(l+@ -yr)x l - %' 

u 1 = -Z + y e 2 -y x 1 + ( py 1 x - x l x i + x e l -xy l + y / e 2 -y/y 2 -x 1 

u 3 = -®e 1 + &x l + ye 3 -yx 3 -l+x l x 2 -% 3 

(8) 

Where %' ,i = 1,2,3 is given as 

Vi 

I 2 I I I 

\X | = A | e 2 | 

(9) 

And A is a diagonal matrix whose diagonals elements diag[C 1 , , ^ 22 , £ 33 ] constitutes the feedback 
coefficients of the controllers, such that 



r* 1 1 [c u o o ir e ,i 

I 2 I I III 

\X 1 = 1 0 ^22 0 || fi 2 | 

U 3 J L« o ^JhJ 



(10) 
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By inserting the elements of (10) in (8) and using (9), the error dynamics reduced to 

(ID 

The Lyapunov stability criteria is used to verify the asymptotic convergence of the error dynamics as t — » oo . 

We choose the Lyapunov function candidate V (.) such thatV (.) > 0; V (.) < 0 . We choose the following 
candidate 

V {e^e ,e ) = ^—{e[ + el + e\) 
20 

(12) 

By using (1 1) in the partial derivative of (12), we have 

V ( ei ,e 2 ,e 3 ) = - £ 21 e\ - £ 

(13) 

Which is Hurwitz for all f > 0 . Consequently, the error dynamics will converge asymptotically to the origin 

as t — > oo 

IV. SIMULATION RESULTS 

The Rikitake dynamo (1), Bullard dynamo (2) and the nonlinear control law (8) were simulated in 
MATLAB environment for the following parameters 0 = 2, ^ = 20,^ = 2,^=4 and 

yr = 2, /r= 3.4641 for initial conditions [x l (0), x 2 (0), x 3 (0)] = [- 1 0, -7 , -4] and 

[^[(O)* y 2 (0), y 3 (0)] = [-16,-5,-1] . The initial conditions of the antisynchronization error dynamics 

becomes[e 1 (0), e 2 (0), e 3 (0)] = [-26, - 12,-5]. The resultant plots are given in the following figures. 




0 0.5 1 1.5 2 0 0.5 1 1.5 2 

t(s) t(s) 

(c) (d) 



www.theijes.com The IJES Page 52 



Anti-Synchronization of the Bullard and Rikitake Dynamo Systems via Nonlinear Active Control 




-10 



0.5 



1 

t(s) 

(e) 



1.5 



Figure 6: Simulated results of the antisynchronized systems - (a) plot of the asymptotically converged error 
dynamics; (b) Plot of the converged control laws; (c) Antisynchronized trajectories xl-yl; (d) 
Antisynchronized trajectories x2-y2 and (e) Antisynchronized trajectories x3-y3. 

V. CONCLUSION 

In this paper, the exponential divergent trajectories of the chaotic Bullard and Rikitake dynamos 
were antisynchronized using nonlinear active control laws. Lyapunov stability criterion was applied to the 
error dynamics inorder to test for asymptotic convergence. The partial derivative of the Lyapunov function 
candidate was Hurwitz, and as a result, the error dynamics and active control laws asymptotically converged in 
transient time. In power system engineering, antisynchronization has found usefulness in mitigation of power 
outage. 
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